Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 7 April 2009 (MN BOgC style file v2.2) 



The Fermi blazars divide 



G. Ghisellini 1 *, L. Maraschi 2 , F. Tavecchio 1 

1 INAF ' — Osservatorio Astronomico di Brera, via E. Bianchi 46, 1—23807 Merate, Italy 
2 INAF - Osservatorio Astronomico di Brera, Via Brera, 28, Milano 1-20121 



ON 

o 
o 

< 



o 

u 

6 
& 

(N 
> 

o 

(N 

cn 
o 
o> 
o 



7 April 2009 



X 



ABSTRACT 

Flat Spectrum Radio Quasars (FSRQs) and BL Lac objects detected in the first three months 
of the Fermi survey neatly separate in the 7-ray spectral index vs 7-ray luminosity plane. 
BL Lac objects are less luminous and have harder spectra than broad line blazars. We suggest 
that this division has its origin in the different accretion regimes of the two classes of objects. 
Using the 7-ray luminosity as a proxy for the observed bolometric one we show that the 
boundary between the two subclasses of blazars can be associated with the threshold between 
the regimes of optically thick accretion disks and of radiatively inefficient accretion flows, 
which lies at an accretion rate of the order of 10~ 2 the Eddington rate. The spectral separation 
in hard (BL Lacs) and soft (FSRQs) objects can then result from the different radiative cooling 
suffered by the relativistic electrons in jets propagating in different ambients. We argue that 
the bulk of the most luminous blazars already detected by Fermi should be characterised by 
large black hole masses, around 10 9 solar masses, and predict that lowering the 7-ray flux 
threshold the region of the a T — L 7 plane corresponding to steep spectral indices and lower 
luminosities will be progressively populated by FSRQs with lower mass black holes, while 
the region of hard spectra and large luminosities will remain forbidden. 
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1 INTRODUCTION 

The Large Area Telescope (LAT) on board the Fermi Gamma Ray 
Space Telescope (Fermi) revealed more than one hundred blazars 
with a significance larger than 10<r in the first three months of op- 
eration (Abdo et al. 2009a, hereafter A09). Of these, 57 are clas- 
sified as flat spectrum radio quasars (FSRQs), 42 as BL Lac ob- 
jects, while for 5 sources the classification is uncertain. Including 

2 radio-galaxies the total number of extragalactic sources amounts 
to 106. Redshifts are known for all FSRQs, for 30 BL Lacs, for 1 
source of uncertain classification and for the two radio-galaxies, 
for a total of 90 objects. Particularly significant is the large num- 
ber of BL Lac objects present in the Fermi sample, which allows 
for the first time an objective comparison of their 7-ray properties 
(spectral shape and luminosity) with those of FSRQs. 

A direct result of the first Fermi blazar sample, discussed by 
A09, is the finding of correlations between their 7-ray and radio 
fluxes/luminosities and between the photon spectral index F 7 in the 
Fermi energy band and their radio luminosity. In the latter plane, as 
noted by A09, blazars seem to define a trend, of increasing T 7 (soft- 
ness) for increasing radio luminosity, with BL Lac objects falling 
at the lower end of this trend. However within each of the two sep- 
arate populations these correlations are not highly significant. 

In this letter we consider BL Lacs and FSRQs as part of a 
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unified blazar population, at variance with the approach of A09. 
Using the 7-ray rather than the radio luminosity (the former be- 
ing more indicative of the bolometric power), we discuss the spec- 
tral index/7-ray luminosity plane where BL Lacs and FSRQs oc- 
cupy different regions. We then propose a simple scenario to ex- 
plain why there is a rather well defined boundary between BL Lacs 
and FSRQs. In our scheme all blazars are essentially similar and 
their different observed properties result from the different accre- 
tion rates feeding their central engines. We show that the observed 
7-ray spectral index dependence on 7-ray (as well as radio) lumi- 
nosity can be naturally explained in the above context. Finally we 
discuss some implications and predictions of the proposed scenario. 
We adopt a cosmology with h = S7a = 0.7 and Qm = 0.3. 



2 THE GAMMA-RAY SLOPE VS LUMINOSITY PLANE 

We derived the 7-ray luminosities for all the blazars with known 
redshift detected by Fermi in the three months survey (excluding 
radio galaxies) using the fluxes named F100 in Table 3 of A09, 
characterising their "average" observed flux (between August 4 and 
October 30 2008). The luminosities were computed according to: 



(i) 
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Figure 1. Energy spectral index vs 7-ray luminosity for all blazars in the list of Abdo et al. (2009a). Blue and red points are BL Lacs and FSRQs, respectively. 
Black symbols correspond to sources already detected in the TeV band (they are all BL Lacs but 3C 279). Labels identify all BL Lac objects, 3C 279 and other 
2 FSRQs. For several BL Lac objects and for the 3 labelled FSRQs we show the observed range of 7-ray luminosity and spectral index, using past EGRET 
or AGILE observations. This is indicated by a solid line segment (dashed when the spectral index is unknown, and we chose the same index as observed by 
Fermi). For 0716+714, the high luminosity point corresponds to AGILE observations (Chen et al. 2008) with a large error on the spectral index (shown by the 
large error bar). The grey stripes at about L 7 = 10 47 erg s — 1 mark the divide between BL Lac objects and FSRQs. 



where is the luminosity distance, a 7 = T 7 — 1, S(yi, Vi) is the 
7-ray energy flux between the frequencies v\ and U2, calculated 
from the photon flux F~,(E > 100 MeV) [ph cm -2 s" 1 ] as: 



1 — a 

S~ 1 (v\,V2) = hv\F^\n{v2/vi) 



Q 7 / 1 



(2) 



Setting v x = 2.42 x 10 22 Hz (0.1 GeV) and u 2 = 2.42 x 10 24 Hz 
(10 GeV), and F 1 = 10" 8 F 7 ,_ 8 ph cm -2 s~\ we have (in units 
of erg cm -2 s _1 ): 



S 7 (0.1,10) 



1.6 x 10 



-12 a 7 -F 7 



[100 1_ 



1 a T 1 



^(0.1,10) = 7.38 x 10" 



a-» = 1 



(3) 



Fig. [T] shows the energy spectral index a 1 given by A09 as 
a function of L 7 between 0.1 and 10 GeV calculated according 
to Eqs. Q}{3] Blue and red symbols correspond to BL Lac objects 
and FSRQs, respectively, while black symbols indicate sources also 
detected in the TeV bancQ The latter are all BL Lac objects except 
for 3C 279, which is an intermediate object with Lya emerging in 
low states (Koratkar, Pian & Urry 1998; Pian et al. 1999). 



1 We have included 0716+714 as a TeV detected BL Lac (Teshima 2008), 
although it is not noted as such in A09. Therefore the total number of Fermi 
bright blazar that have been also detected in the TeV band is 12. 
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Blazars are highly variable sources, especially at high ener- 
gies. To illustrate the possible range of variability, for a few sources 
we connect with a line segment the spectral index and 7-ray lumi- 
nosity observed by Fermi to the values observed in the past by the 
EGRET instrument onboard the Compton Gamma Ray Observa- 
tory satellite (Nandikotkur et al. 2007), or observed by the AGILE 
satellite. Note that variable sources often "move" orthogonally to 
the correlation defined by the ensemble of sources, i.e. they become 
harder when brighter (with the exception of 3C 454.3). 

Moreover, for some of the sources detected in the TeV range, 
we can infer the level of 7-ray activity during high states of TeV 
emission (assuming that the flux in the Fermi band would have been 
about the same as the TeV flux). For these sources we indicate the 
range of luminosity variability by dotted line segments, maintain- 
ing the spectral index measured by Fermi for both states. 

As can be seen, the high and the low 7-ray states of single 
sources can be dramatically different. This implies that the distri- 
bution in luminosity within each blazar class is largely affected by 
the variability of the sources. 

The a 7 -L 7 plane (Fig.[]} strongly suggests a separation of 
the two groups in regions defined by a 7 > 1.2, L 1 > 10 47 erg s _1 
for FSRQ and a-, < 1.2, L-, < 10 47 erg s _1 for BL Lacs. At the 
same time there is continuity between the two blazar subclasses in 
that their properties overlap at intermediate values of both spectral 
index and 7-ray luminosity. 

The separation in spectral indices, already evident in Fig. 9 of 
A09, occurs at a 7 ~ 1.2 Only 5 FSRQs have spectra flatter than 
1.2, and only 5 BL Lacs have spectra steeper than 1.2. The LAT 7- 
ray data have been systematically analysed by A09 fitting a single 
power law model to the spectral data of each source. In many cases 
the intrinsic spectral shape could be more complex (i.e. a curved 
spectrum). Fig. 10 of A09 shows as illustrative examples the 7-ray 
SEDs of 3C 454.3 (energy index a 7 = 1.41±0.02), AO 0235+164 
(a 7 = 1.05 ± 0.02) and Mkn 501 (a 7 = 0.70 ± 0.09). In the first 
two cases the data are not well represented by the single power- 
law model, showing significant curvature and indicating SED peaks 
below 100 MeV for 3C 454.3 and between 1 and 10 GeV for AO 
0235+164, while Mkn 501 shows a very hard 7-ray SED implying 
that the SED peak should be above 100 GeV. It is therefore reason- 
able to assume that the spectral index measured by Fermi is a proxy 
for the peak energy of the 7-ray SED, at least in a statistical sense. 
Thus the 7-ray SEDs of most FSRQs, with a 7 > 1.2, should peak 
at or below the low energy range of the Fermi energy window, at 
lower peak energies than those characterising most of BL Lacs ob- 
jects. As expected, BL Lacs detected in the TeV band tend to be 
harder than other BL Lacs (with the exception of BL Lac itself), 
indicating a high peak energy, in turn necessary to detect them in 
the TeV band. 

The separation in 7-ray luminosity between BL Lac objects 
and FSRQs is also striking. There are only 3 BL Lac objects (PKS 
0537-441 with z = 0.892, AO 0235+164 with z = 0.94 and PKS 
0426-380 with z = 1.12) with L 7 > 10 47 erg s" 1 . The redshifts 
of these 3 sources are among the highest of the BL Lac sub-sample 
(see Fig. 16 of A09). Moreover, all three sources do have broad 
lines, visible in their low emission states (see Sbarufatti et al. 2005 
for PKS 0426-380; Pian et al. 2002 for PKS 0537-441; Raiteri et 
al. 2007 for AO 0235+164). 

In principle, the two blazar classes could be separated because 
of a different average Doppler boosting, but all indications so far 
for the 7-ray emitting BL Lacs and FSRQs suggest instead similar 
values of the Doppler boosting, with the exception of TeV emitting 
BL Lacs, requiring a stronger boost (and thus a smaller viewing 



angle). These inferences come from superluminal motion (Jorstad 
et al. 2001; Kellerman et al. 2004), direct modelling of their SED 
(e.g. Konopelko et al. 2003; Finke, Dermer & Bottcher 2008) and 
very rapid TeV variability (e.g. Begelman, Fabian & Rees 2008; 
Ghisellini & Tavecchio 2008b). Thus the idea that the luminosity 
sequence of BL Lacs and FSRQs is due to larger viewing angles 
for BL Lacs is disfavoured. 

2.1 The blazar sequence in 7-rays 

The SEDs of blazars show two main components: the first one, 
peaking between the far infrared and the X-ray band, is unani- 
mously attributed to synchrotron emission; the second, peaking be- 
tween 10 MeV-100 GeV, is widely attributed to the Inverse Comp- 
ton process involving the same electrons producing the synchrotron 
component (e.g. Maraschi et al. 1992; Bloom & Marscher et al. 
1996; Sikora et al. 1994; Dermer & Schlickeiser 1993). 

Fossati et al. (1998) and Ghisellini et al. (1998) proposed 
that the SEDs of blazars form a "spectral sequence" whereby the 
peaks of the two components are governed by luminosity, and are 
at higher energies for lower luminosities^. This was interpreted in 
terms of a stronger radiative cooling for more powerful sources, re- 
sulting in smaller energies of the electrons emitting at the peaks. As 
discussed above, a high energy peak above ~10 GeV corresponds 
to a hard spectrum in the Fermi range (rising, in vF v ), while a 
peak frequency below ~1 GeV corresponds to a soft spectrum in 
the Fermi range. Therefore the trend apparent in the a 7 -L 7 plane 
strongly supports the "sequence" concept, which was proposed on 
the basis of radio and X-ray selected samples, where only a frac- 
tion of objects had 7-ray data. Fig.[T]represents the 7-ray selected 
version of the blazar sequence. 

Single objects, instead, often behave in the opposite way (i.e. 
harder when brighter). This can be explained, if, in single objects, 
an unchanged cooling faces a fast increasing heating (i.e. injection 
of more energetic electrons), shifting the typical electron energies 
to larger values. In the next section we present a possible interpre- 
tation of the separation of FSRQs and BL Lacs in the a 7 -L 7 plane. 



3 INTERPRETATION OF THE BL LAC - FSRQ DIVIDE 

The neat separation between BL Lacs and FSRQs in the a 7 -L 7 
plane together with their close contiguity are really remarkable. 
A possible interpretation could simply be to assume two indepen- 
dent populations with different physical conditions as to the ori- 
gin, propagation and radiative properties of their jets. This however 
would not account for the contiguity which should be coincidental. 
Figs. 1 1 and 17 of A09 (a 7 vs z and L 7 vs z) indicate that the sep- 
aration and the contiguity behaviour of the two blazar sub-classes 
exist also in these planes. 

Building on previous work on the systematic properties of 
blazars we suggest that this division has a physical origin, related to 
the different mass accretion rates onto the central engine. We first 
summarise some preliminary steps as follows: 

(i) Radiative models allow to estimate the total jet power Pj car- 
ried in the form of bulk and internal energy of protons, electrons 
and magnetic field. For powerful FSRQs a fraction e ~ 0.1 of Pj 

2 The original idea considered the luminosity as the main parameter gov- 
erning the blazar sequence, while the refined version (Ghisellini & Tavec- 
chio 2008a) considers the luminosity in Eddington units. 
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is radiated at short distances (i.e. hundreds of Schwarzschild radii) 
from the central black hole. This fraction can be larger for BL Lacs 
(Celotti & Ghisellini 2008). 

(ii) In powerful FSRQs, Pj is about one order of magnitude 
larger than the observed accretion disk luminosity Ld (Maraschi 
2001; Maraschi & Tavecchio 2003; Celotti & Ghisellini 2008; 
Ghisellini & Tavecchio 2008a; see also the earlier results of Rawl- 
ings & Saunders 1991). For a "standard" Shakura & Sunyaev 
(1973) accretion disk, Ld = r]Mc 2 , where M is the accretion rate, 
with rj ~ 0.1. It follows that the jet power Pj is of the same order 
as the accretion power Mc 2 . This is also supported theoretically (at 
least at the "dimensional level") since, even if the black hole rota- 
tion played an essential role as in the Blandford & Znajek (1977) 
model, accretion would still be needed to provide the magnetic field 
necessary to extract the rotational energy. We then assume that the 
(order of magnitude) equality Pj ~ Mc 2 has general validity for 
strong jets, being aware of the significant uncertainties associated 
with it. 

(iii) Finally, the 7-ray luminosity dominates the electromag- 
netic output of powerful FSRQs and is at least a sizeable fraction of 
it in less powerful BL Lacs. The observed L 7 ~ T 2 ePj (the factor 
F 2 accounts for light aberration). 

Adopting the relations above, it then follows that the max- 
imum observed L 7 should correspond to the largest black hole 
masses accreting at the maximum, Eddington, rate 

L 7 ,max ~ er 2 M Edd c 2 ~ 1.5 x 10 47 er 2 Af 9 erg s" 1 (4) 

where Mg is the mass of the black hole in units of 10 9 solar masses. 
The 7-ray luminosity of the most luminous FSRQs (Fig.[TJ is of the 
order of several xlO 48 erg s _1 , corresponding to eF 2 M9 ~ 30. 

Fig. [T] shows a divide at a 7-ray luminosity roughly a factor 
100 below the maximum (grey stripes). Maintaining the same value 
of er 2 M 9 ~ 30, we conclude that the jet power, hence the accre- 
tion rate, is a factor ~100 below the maximum. Therefore, at the 
dividing luminosity, the accretion rate is of order: 

Mdiv « 10- 2 M Edd (5) 

Crucial, in this simple derivation, is the assumption that the most 
luminous FSRQs and the most luminous BL Lacs (excluding the 
three "outliers" mentioned above) have the same factor eT 2 M, and 
in particular the same black hole mass M. This is likely, since we 
are dealing with the most luminous objects of the two classes (this 
also applies to the original blazar sequence). Indeed, for TeV BL 
Lacs, which span the entire region of BL Lacs, the available es- 
timates indicate high black hole masses (e.g. Wagner 2008, Ghis- 
ellini & Tavecchio 2008a), in some cases exceeding 10 9 i\/©. In the 
case of FSRQs (more luminous than BL Lacs) Fermi could have 
already started to explore smaller black hole masses. In fact, the 
object PMN 0948+0022 (a 7 = 1.6 ± 0.14 and L 7 ~ 10 47 erg 
s _1 ) is a NLSyl, with a black hole mass around 1O 8 M0 (see Abdo 
et al. 2009b and Zhou et al. 2006). This easily accounts for the few 
FSRQs with L 7 < 10 47 erg s _1 (i.e. they should have black hole 
masses smaller than the rest of FSRQs). 

The derived value of Mdiv is naturally interpreted as due to the 
transition between different accretion regimes: for high accretion 
rates, in Eddington units, a standard optically thick geometrically 
thin, radiatively efficient Shakura-Sunyaev (1973) disk is formed 
producing a luminous Broad Line Region through which the jet 
must propagate suffering strong radiative losses. For M ~ M-^dd 
some of the disk luminosity could be advected into the black hole, 



resulting in r\ < 0.1, but this would not affect the presence of the 
BLR, nor the power of the jet. 

When M <C M^dd the accretion flow becomes optically thin 
and radiatively inefficient and could be described as an ion sup- 
ported torus (Rees et al. 1982) or an advection dominated accretion 
flow (ADAF, see e.g. Narayan, Garcia & McClintock 1997), an adi- 
abatic inflow-outflow (ADIOS, Blandford & Begelman 1999) or 
a convection dominated flow (CDAF, Narayan, Igumenshchev & 
Abramowicz 2000). In the latter case despite the uncertainties in 
the dynamics of the flow the surrounding region will certainly be 
deprived of photons with respect to the high accretion rate regime, 
thus justifying the absence of broad emission lines and the reduced 
radiative losses for the propagating jet. 

Since FSRQs are associated with radio sources with FR II 
morphology while BL Lac objects are usually found in FR Is, the 
same critical rate should separate the two parent populations. And 
indeed the finding by Ledlow & Owen (1996) that the radio lu- 
minosity threshold between the two classes depends on the optical 
magnitude of the bulge of the host galaxy (i.e. on the mass of the 
central black hole), was interpreted by Ghisellini & Celotti (2001) 
in terms of a critical accretion rate, coincident in value with what 
found here (see also more recent work by Xu, Cao & Wu 2009). 



4 DISCUSSION AND CONCLUSIONS 

The spectral shape/7-ray luminosity plane revealed a clear sepa- 
ration between FSRQs and BL Lac objects. Moreover, there is a 
well defined trend, with BL Lac objects being harder and less lumi- 
nous than FSRQs. At the first order, this strongly confirms the idea 
that blazars obey a spectral sequence, controlled by the bolometric 
luminosity, of which the 7-ray one is good proxy. 

The other striking feature is the little overlap, in 7-ray lumi- 
nosity, of the two classes of blazars. We suggest that this has a phys- 
ical origin, and have interpreted it in terms of a critical accretion 
rate Mdiv Below Mdiv accretion becomes radiatively inefficient, 
the corresponding ionising flux becomes weaker, making weaker 
(or absent) broad emission lines. These are BL Lac objects. Their 
jets thus propagate in a medium starved of external radiation (weak 
disk, weak lines), and this makes the emitting electrons accelerated 
in the jet to cool less, and to reach very high energies. Their emitted 
high energy spectrum, produced mainly by the synchrotron self- 
Compton mechanism, is less luminous and harder than in FSRQs. 
The latter sources in fact have a radiatively efficient accretion disk 
and a corresponding "standard" broad line region. If the jet dissi- 
pates most of its power within the broad line region, then the emit- 
ting electrons will efficiently cool (mostly by external Compton), 
will reach only moderate energies, and will produce a high energy 
peak below 100 MeV. This scenario is in perfect agreement with 
the analogous division proposed for radio-galaxies, that in turn are 
though to be the progenitors of BL Lacs (FR I) and FSRQs (FR II). 

According to this idea the power of the jets in all blazars is 
proportional to the accretion rate, that can change in a continuous 
way from one objects to the other (and even in single objects). The 
observed discontinuity between BL Lacs and FSRQs is the result 
of a discontinuity in the accretion regime occurring within a rather 
narrow range of M (from standard to radiatively inefficient). Thus, 
despite the different "look" of FSRQs and BL Lacs, they all belong 
to a sequence of jet powers. From this point of view, it is therefore 
reasonable to treat them as a single class. 

Evidences that Pj is of the same order of Mc 2 can be traced 
back to Rawlings & Saunders (1991), and since then confirmed by 
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several works (e.g. Celotti et al. 1997; Cavaliere & D'Elia 2002; 
D'Elia, Padovani & Landt 2003; Maraschi & Tavecchio 2003; Sam- 
bruna et al. 2006; Allen et al. 2006; Celotti & Ghisellini 2008; 
Kataoka et al. 2008). Also, the discovery that large scale jets of 
powerful blazars are strong X-ray emitters, detectable and resolv- 
able by the Chandra satellite, offered a novel way to measure Pj , 
confirming that in these sources the jet power is not only of the 
same order than Ale 2 , but it is also constant all along the jet (Tavec- 
chio et al. 2004; Cheung et al. 2006). 

Also the idea that the differences between powerful FSRQs 
(and their FR II radio-galaxy counterpart) and BL Lacs (and FR 
I radio-galaxies) is due to a difference in the accretion regime of 
their disks has been already put forward (Ghisellini & Celotti 2001 ; 
Jester 2005) and used to explain the different cosmic evolution of 
these two classes of blazars (see e.g. Bottcher & Dermer 2002; Cav- 
aliere & D'Elia 2002). Interestingly, also in radio-quiet AGNs there 
is evidence that the relative strength of the broad emission lines 
is a function of the accretion disk luminosity in Eddington units, 
with the possibility of a very weak BLR for accretion rates below 
10 -2 Z/Edd (for review, see Ho 2008). What is new, in this paper, 
is the interpretation, in terms of accretion rates, of the intriguing 
separation of FSRQs and BL Lacs detected by Fermi in the a~ l -L 1 
plane. A (likely) cosmic evolution of A-/, while not affecting our 
results, would account for the existing hints of a negative and pos- 
itive cosmic evolution of BL Lacs and FSRQs, respectively (see, 
e.g. A09). 

There are a few checks to put our scenario on test. Perhaps the 
simplest is to have clues on the black hole masses of the most lumi- 
nous BL Lacs and FSRQs. We predict them to be of the same order. 
A related check is to directly find the accretion disk luminosity for 
the most luminous FSRQs. If also the black hole mass is known or 
derived, then we can test if their disk are indeed emitting close to 
the Eddington limit, as we have assumed. 

Although exceptions exist for any rule, it should be interesting 
to study in more detail the three BL Lacs that are "outliers", lying 
in the FSRQ-region. Since they do have broad lines, we propose 
that they are normal FSRQs whose non-thermal continuum is so 
enhanced by beaming that it hides the broad lines. If this is the case, 
the FSRQs/BL Lacs division, now in terms of the line equivalent 
width, could be substituted by a more physical division in terms of 
broad line luminosity. 

In just 3 months of all sky survey (and selecting only the 
brightest ones), Fermi has already doubled the number of known 
7-ray emitting blazars. Soon we will have to deal with thousands 
of them. Since blazars are very variable, perhaps spending a limited 
amount of time in high state, continuing to patrol the all sky will 
surely reveal exceptional states of sources otherwise unnoticeable. 
On the other hand, the increase in sensitivity due to the increased 
exposure time will allow to find many less powerful jets and less lu- 
minous 7-ray luminosities. The divide will then be not as clear as it 
is now, since many FSRQs (with black hole masses below 1O 9 M ) 
will occupy the region below L 1 — 10 47 erg s -1 . On the other 
hand, the prediction is that no or very few BL Lacs will be found 
above this limit. These improved samples will surely be useful to 
confirm, or reject, or constrain our scenario. 
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